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To understand the link between doping and electronic 
properties in high temperature superconductors, we report 
first-principles calculations on the oxygen doping effect for 
the single layer cuprate HgBa2Cu04+i. The doping effect is 
studied both by supercell and single cell virtual crystal type 
approaches. We find ionic behavior of the dopant atom up 
to an optimal doping concentration oi S — 0.22. The excess 
oxygen attracts electrons from the Cu02 plane leading to an 
increase of the hole concentration in this building block. The 
maximum amount of holes is reached when the oxygen p-shell 
is nearly completely filled. As a consequence the density of 
states at the Fermi level exhibits a pronounced maximum at 
the optimal oxygen content. The calculated hole content as 
a function of doping is in accordance with experimental find- 
ings. 



It is commonly believed that the charge carriers re- 
sponsible for the superconducting pairing of high Tc 
cuprates are holes mainly confined to the Cu02 layers. 
Similar correlations between Tc and and ng/m*, with Hs 
and TO* being the charge carrier density and the effective 
mass, respectively, do not only include high Tc materials, 
but also other families of superconductors A univer- 
sal relationship between Tc/T^"-^ and the hole content p, 
exhibiting a common feature among the p-type high-Tc 
superconductors, was presented by Zhang et al. Q. This 
dependence is characterized by a plateau, where the vari- 
ation of Tf./T^"'^ is independent of the compound con- 
sidered. 

In all these materials of interest the actual amount 
of holes is driven by doping. However, for most of the 
theoretical models describing the superconducting phase 
transition the hole content rather than the doping level is 
the crucial input parameter. Above all, hardly anything 
is known about the relationship between these two physi- 
cal quantities. Without a detailed knowledge of how dop- 
ing influences the number of carriers in the normal state, 
also the superconducting properties will lack a profound 
understanding. In view of being able to tune the amount 
of carriers and thereby the superconducting transition 
temperature a full clarification of how doping affects the 
electronic structure is not only highly interesting, but 
even inevitable. In this context the following questions 
are raised: (i) Where does the ecxess charge go upon 
doping? (ii) How does doping influence the carrier con- 
centration in the copper-oxygen planes? (iii) What limits 



the amount of holes in these building blocks? (iv) How 
does the density of states (DOS) behave as a function of 
doping? 

In order to address these topics, we have carried out 
electronic structure calculations based on density func- 
tional theory (DFT). For that purpose we have chosen 
the simplest representative of the Hg-based high Tc com- 
pounds, the Cu02 single-layer HgBa2Cu04+5 (Hgl201). 

This family of layered structures with the formula unit 
of HgBa2Ca„_iCu„02„+2+5 (H represents the class of su- 
perconducting materials with the highest transition tem- 
perature so far. It reaches a value of 136 K for n — 2> 
which can be even enhanced by applying pressure [Q. 
Upon adding some excess oxygen, Tc can be dramati- 
cally increased For Hgl201, a systematic study over 
a wider range of doping (0.03 < 5 < 0.4) ||^ revealed 
a parabolic-like Tc dependence on the oxygen content 5 
exhibiting a maximum at 5opt — 0.22. At this optimal 
doping concentration Tc is 97 K for ambient pressure, 
and climbs up to 118 K at 24 GPa, which is the high- 
est Tc of any single-layer cuprate. Although the mercury 
based high Tc compounds are not explicitely included in 
the experimental analysis described above ||l|,|| one can 
expect them to fit into this framework. 

Theoretically, doping effects in this class of materials 
have been studied for a high doping level of (5 = 1/2 
only [^,0 revealing strong covalent bonding in the basal 
plane with large overlap of the Hg and the dopant oxygen 
orbitals. This oxygen concentration, however, has never 
been achieved experimentally. Therefore a detailed study 
of the much more important low doping regime is highly 
desirable. In our investigations special emphasis is put 
to the doping induced redistribution of the hole concen- 
tration in the Cu02 layer. In addition, the doping effect 
on the crystal structure is studied. In order to make ar- 
bitrary doping levels feasible we treat an excess oxygen 
concentration 5 by adding the appropriate amount of va- 
lence electrons 6 x (5 (6 electrons for the fully doped cell) 
to the crystal in a single-cell calculation. To guaran- 
tee charge neutrality the corresponding positive charge 
was placed in a sphere at the doping oxygen site, i.e. 
in the middle of the Hg squares in the basal plane. As 
we will see below, this virtual crystal approach is jus- 
tified by the ionic behavior of the dopant oxygen. We 
have cross-checked our results by supercell calculations 
for doping levels of <5 = 1/6 and 5 = 1/4, and in addition 
for 5 = 1/2. For the systematic study of the doping- 
dependent electronic structure, the crystalline data as a 
function of doping were taken from Ref. S. In a next 
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step, we also investigate the effect of doping on the lat- 
tice parameters and atomic positions by total-energy and 
atomic-force calculations. 

All calculations are carried out within the full-potential 
linearized augmented plane-wave (LAPW) method ||] 
utilizing the WIEN97 code |0. Exchange and corre- 
lation effects are treated within the local density approx- 
imation (LDA). The Brillouin zone integrations are car- 
ried out on a 20 X 20 X 8 k-point mesh (consisting of 
220 points in the irreducible wedge) applying a Gaussian 
smearing of 0.002 mRy. Our basis set includes approxi- 
mately 1500 LAPW's supplemented by local orbitals for 
the low- lying semicore-states Hg-5p, Gn-Sp and 0-2s. 
We used atomic sphere radii of 2.0, 2.2, 1.9, and 1.55 a.u. 
for Hg, Ba, Cu and O, respectively. The atomic- like ba- 
sis functions used within the LAPW method allow an 
analysis and the orbital symmetry decomposition of the 
electronic charge within the atomic spheres. For the la- 
belling of the oxygen atoms see Fig. ??. 

As already found by previous calculations carried out 
for the undoped material [^,0, the Fermi level inter- 
sects a single free-electron-like two-dimensional half-filled 
dpa* band, with its states being of Cu(dj,2_j^2) and 0{px) 
character. Upon doping the effect on the electronic 
structure is twofold. First, the dopant adds states at 
the Fermi level, where these new charge carriers cause 
a shift of Ep. Second, the shift of the Fermi level 
changes the carrier concentration in the Cu02 plane. 
This can be seen in Fig. |^ where the corresponding 
losses in the partial charges with respect to the undoped 
system are displayed. The occupation number for the 
copper dx2-y2 orbital Q(dx2_y2) is decreased by 0.09 e 
(from 1.443 e for the undoped case to 1.353 e for 5opt), 
and similarly the Ol(px) charge Q{px) drops by 0.030 
e (from 0.982 to 0.952). Both orbitals clearly exhibit 
their minimum occupation number (maximum amount 
of holes) as a plateau around the experimentally ob- 
served optimal doping content. The hole concentration 
within the copper-oxygen planes with respect to the un- 
doped case is obtained by [Q{dx2_y2){Q) + 2 x Q{px){^)] - 
[Q{dx2-y2){6) + 2x Q{px){S)] also displayed in Fig. |l|. In 
excellent agreement with experiment |^ its optimal value 
is found to be 0.16 e. 

The increase of the hole content can be generally un- 
derstood in an ionic picture: The excess oxygen attracts 
electrons from the copper-oxygen plane thereby increas- 
ing the hole concentration in this region. In fact, the 03 
states do not show strong overlap with the neighboring 
Hg states in the entire doping regime up to optimal dop- 
ing. This picture is confirmed by supercell calculations 
for doping concentrations of (5 = 1/6 and 5=1/4 (with 
one excess oxygen 03 per six and four single cells, respec- 
tively). It shows a nearly flat oxygen band at Ep with 
vanishing Hg-d admixture, whereas the Hg-d states are 
found somewhat below. The total and partial densities 
of states for (5=1/4 are shown in Figure It exhibits 
covalent bonding of Cu and 01 with peaks at -0.35 eV 
and a small feature right below Ep. Covalent bonds are 



also formed by Hg and the apical oxygen 02, whereas the 
dopant atom 03 has the highest contribution to the DOS 
at the Fermi level with hardly any common feature with 
mercury. This fact is supported by the partial charges 
in the Hg sphere which are not affected by doping. Only 
in the overdoped region, investigated for 5 = 1, and for 
S = 1/2 by a supercell calculation, stronger bonds be- 
tween oxygen and mercury atoms are formed as was also 
found in the three-layer compound for the heavily oxy- 
genated material 0. 

Let us turn now to the question what the limiting fact 
for the creation of holes in the Cu02 plane is. We pro- 
pose that the hole concentration reaches its maximum 
when the dopant band is nearly completely occupied, i.e. 
for a closed 03 shell. Since the excess oxygen attracts 
electrons from the Cu02 plane, its orbitals will be sys- 
tematically filled up not only by the electrons provided 
by the oxygen atom itself, but also by those coming from 
the Cu02 unit. In order to quantitatively prove our hy- 
pothesis, we have repeated our calculations with a muffin 
tin radius of 2.5. a.u. for the artificial 03 atoms with 
nuclear charge oi Z ~ 6(5. This sphere is big enough to 
include the doping charge and therefore allows the inter- 
pretation of the charges inside this sphere as absolute oc- 
cupation numbers. We indeed find a nearly closed shell 
at the experimentally observed optimal doping concen- 
traion Sopt — 0.22 which limits the amount of holes 
in the Cu02 plane an therefore most probably the su- 
perconducting transition temperature. Furthermore, the 
occupation of the 2pz shell of the doping O atom (full 
oxygen atom here) is found by 4cell calculation. We used 
a muffin tin sphere pf 3.0 a.u. in this calculation and an O 
anion is found with the negative charge of 1.4 e. Further 
increase of the muffin tin sphere was not possible due to 
computational difhculties. Therefore, at this moment it 
is not easy to estimate whether a anion is formed 
or an anion with less then 2.0 e negative charge. Never- 
theless, the presence of anion dopant state at the doping 
site can be clearly seen on the basis of LDA calculations. 
In the overdoped region covalent bonds between 03 and 
Hg are formed. By this change the pure ionic charac- 
ter of the dopant oxygen gets lost and therefore no more 
carriers are attracted from the Cu02 plane which results 
in the broad plateau on Fig. in accordance with the 
experimental findings |^]. 

These findings also show up in the density of states 
which in the optimal doping regime exhibits an almost 
filled band right below Ep. The calculated DOS at the 
Fermi level N{Ep) as displayed in Fig. |^ has a pro- 
nounced maximum at 5q = 0.22. It contains large con- 
tributions from the Cu02 plane, but also a considerable 
amount of apical oxygen and some 03. All the contri- 
butions peak at the optimal doping content. ^From the 
density of states at Ep the electronic specific heat coeffi- 
cient can be determined. Assuming a similar behavior for 
other high- Tc materials like the Bi-based compounds, the 
measured specific heat coefficient 7 of Bi2Sr2CaCu208+5 
as a function of doping [|l^ can be understood by the 
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pronounced increase of the DOS around optimal dop- 
ing. Moreover, the calculated maximum value of 7=4.4 
mJ/(mol K^) is in excellent accordance with experimen- 
tal values for other high Tc materials at optimal doping 

As already mentioned, we crosschecked our results with 
supercell calculations for 5 = 1/6 and 5 — 1/4 in order 
to ensure that no artifacts may arise due to our single- 
cell approach for the treatment of doping. The partial 
charges in all the atomic spheres agree very well with 
the corresponding single-cell calculation (within one hun- 
dredth of an e). E.g. for 5 — 1/4 the copper Q{dx^-y2) 
partial charge is 1.364 compared to 1.353 in the single- 
cell approach, the oxygen Q{px) charges are 0.949 and 
0.952, respectively. For 6=1/6 the deviations are also 
below 1%. Excellent agreement is also found for the den- 
sity of states at the Fermi level, as can be seen for both 
cases in Fig. ^. 

Not only the electronic properties, but also the crys- 
tal structure is influenced by doping. The redistribution 
of bond lengths can again be explained within the ionic 
picture. Since doping does not change the symmetry of 
the crystal, it only leads to displacements of Ba and the 
apical oxygen in z-direction. With increasing oxygen con- 
tent at the doping site, the positive Ba ions are expected 
to be attracted by the oxygen, whereas the negatively 
charged 02 ions should be repelled. Indeed a shift of the 
Ba position towards the basal plane is observed together 
with a displacement of the 02 ions away from this plane. 
In Fig. ^ the corresponding z-coordinates of both atoms 
are shown which were calculated by fully optimizing the 
crystalline data for selected doping concentrations: For 
6=0, 0.08, 0.17, and 0.23 the unit cell volumes, c/a ra- 
tios, and atomic positions were obtained by total-energy 
and atomic-force calculations using a similar procedure 
like in Ref . (m) . The results are displayed in Table |. All 
trends are in agreement with experimental observations 

: These are a shrinkage of the unit cell upon doping, a 
shortening of the lattice parameter a, an increase of c for 
small doping levels followed by a nearly linear decrease, 
an increasing distance between the basal plane and the 
Ba layer, as well as a reduction of the distance between 
the layers of Ba and 02. In agreement with measured 
data, the latter two quantities show a less pronounced 
dependence up to (5 = 0.08 and a stronger one for higher 
oxygen contents. 

In summary, we have investigated the relationship be- 
tween doping and electronic and structural properties 
in the high temperature superconductor IIgBa2Cu04+5. 
We find that for doping concentrations up to optimum 
the dopant does not form covalent bonds with the neigh- 
boring Hg atoms, but shows purely ionic behavior. The 
doping oxygen adds states close to the Fermi level show- 
ing nearly pure 03 character. Due to its high electronega- 
tivity the excess oxygen attracts electrons from the Cu02 
plane resulting in a redistribution of charge carriers and 
an increase of the number of holes in this plane. At 
optimal doping the DOS at the Fermi level exhibits a 



pronounced maximum in accordance with experimental 
specific heat data of other cuprates. The doping-induced 
charge redistribution is accompanied by a shrinkage of 
the unit cell and displacements of Ba and the apical oxy- 
gen 02 towards and from the basal plane, respectively. 
The amount of holes created in the Cu02 plane is limited 
by the fact that at optimal doping a closed oxygen shell 
is reached. We conclude that this behavior is a general 
feature among different families of high Tc compounds. 
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TABLE I. Crystalline data as a function of doping ob- 
tained from total-energy and atomic-force calculations. 
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a [A] 


c[A] 


[A3] 


z{Ba) [c] 


zi02)[c] 


0.00 


3.8574 


9.7228 


144.6676 


0.3003 


0.2014 


0.08 


3.8192 


9.8191 


143.2211 


0.2965 


0.2023 


0.17 


3.8058 


9.6887 


140.3283 


0.2851 


0.2055 


0.23 


3.8312 


9.5604 


140.3280 


0.2808 


0.2074 



FIG. 3. Total density of states at Ef per eV and formula 

unit for HgBa2Cu04+f as a function of doping. The right axis 
indicates the electronic linear specific-heat coefficient 7 in mj 
per mol and K^. The maximum value is found exactly at the 
experimentally observed optimum doping concentration. The 
results for the supercell calculations are indicated by stars. 



FIG. 1. Doping induced changes in the partial charges 

within the atomic spheres of Cu (half filled circles) and 01 
(open circles) and total amount of holes (filled circles) in the 
Cu02 plane as a function of the doping concentration S. 



FIG. 4. Doping induced changes in the atomic positions of 
Ba and 02 in units of the lattice parameter c as function of 
doping. 



FIG. 2. Total and site-projected densities of states in states 
per eV and unit cell volume Q for (5 = 0.25 obtained by a 
supercell calculation containing one 03 atom per four formula 
units. (The Fermi level is set to zero.) 
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